jordan_square <-
function(x){
	out <- (lambda1(x)^2)*u1(x) + (lambda2(x)^2)*u2(x)
	return(out)
}

